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We demonstrate that the current-voltage (I-V) characteristics of resistively and capacitively 
shunted Josephson junctions (RCSJ) hosting localized subgap Majorana states provides a phase 
sensitive method for their detection. The TV characteristics of such RCSJs, in contrast to their 
resistively shunted counterparts, exhibit subharmonic odd Shapiro steps. These steps, owing to their 
subharmonic nature, exhibit qualitatively different properties compared to harmonic odd steps of 
conventional junctions. In addition, the RCSJs hosting Majorana bound states also display an 
additional sequence of steps in the devil staircase structure seen in their I-V characteristics; such 
sequence of steps make their I-V characteristics qualitatively distinct from that of their conventional 
counterparts. A similar study for RCSJs with graphene superconducting junctions hosting Dirac- 
like quasiparticles reveals that the Shapiro step width in their I-V curves bears a signature of the 
transmission resonance phenomenon of their underlying Dirac quasiparticles; consequently, these 
step widths exhibit a tt periodic oscillatory behavior with variation of the junction barrier potential. 

We discuss experiments which can test our theory. 

PACS numbers: 05.70.Ln, 05.30.Rt, 71.10.Pm 


I. INTRODUCTION 

The possibility of realization of Majorana zero inodes, 
particles with anyonic statistics described by real wave- 
functions, has attracted tremendous interest in recent 
yearjiJ. Several suggestions regarding condensed matter 
systems which can host such fermions have recently been 
put fortlPHi. Out of these, the most promising ones for 
experimental realization turns out to be those which host 
Majorana modes as localized subgap states in their super¬ 
conducting ground stat^^^^^. Typically, the occurrence of 
such states require unconventional superc onduc ting pair¬ 
ing symmetry such as p- or d-wave pairin^i^lUJ However 
recent proposals have circumvented this requirement; it 
was shown that such bound states can occur either at the 
end of one-dimensional (ID) wire in a magnetic held with 
spin-orbit coupling an d in the presence of a proximate 
s-wave superconductoJ^ or in superconducting junc¬ 
tions atop a topological insulator surface hosting Dirac 
fermions on the surface of a topological insulators^. Such 
Majorana fermions leave their signature as a midgap peak 
in tunneling conductance measuremenlP^ as well through 
fractional Josephson effeclP^. 

Another interesting phenomenon in recent years has 
been the discovery of materials whose low-energy quasi¬ 
particles obey Dirac-like equations. These materials 
are commonly dubbed as Dirac materials; graphene 
and top ologica l insulators are common examples of such 
material j^'^bsi These materials can exhibit superconduc¬ 
tivity via proximity effect with Cooper pairing o ccurri ng 
between Dirac electrons with opposite momentun J^^b7 t 
is well known that transport properties of such supercon¬ 
ductors differ from their conventional counterparts and 
can serve as experimental signat ures o f the Dirac nature 
of their constituent quasiparticleji^E^. 

The experimental detection of Majorana modes have 


mainly relied on measurement of either midgap peak^ 
or detection of even Shapiro steps in a Josephson junc¬ 
tions of superconductors hosting Majorana modej^^. The 
effectiveness of the former set of experiments in detec¬ 
tion of Majorana modes has been questioned since the 
midgap peak did not lead to the expected 2e^//i value 
of the tunneling conductance and could have also oc¬ 
cur from several other effects such as presence of mag¬ 
netic impurities leading to Kondo effectP^ and impurity 
induced subgap states^. In this sense, the presence of 
even Shapiro steps at U = nhwj/e (and the absence of 
odd ones at (2n -|- l)Hijjj/2e) in Josephson current mea¬ 
surement, where ujj is the Josephson frequency, n is an 
integer, and V is the applied external voltage, provide a 
more dehnite detection of such fermions since they con¬ 
stitute a phase sensitive signature which is free of ef¬ 
fects of disordeili^. Consequently, theoretical studies of 
AC Josephson effect for unconventional superconductors 
which hos ts Ma jorana modes has received a lot of atten¬ 
tion lateljpmni Theoretical studies of Josephson effect 
in graph ene Josephson junctions has also been carried 
oufJUMJ T was shown that the critical current of such 
junctions shows a novel oscillatory dependence on the 
barrier potential of the junctions. However, the features 
of AC Josephson effects in either of these systems for a 
resistively and capacitively shunted Josephson junction 
(RCSJ) in the presence of an external radiation has not 
been studied previously. In this context, we mention 
that the analysis presented in the paper for junctions 
hosting Majorana subgap states is somewhat idealized 
in the sense that it does not provide a full treatment 
of Landau-Zenner tunneling and other related quantum 
effects; however we do provide a qualitative discussion 
and identify a regime of junction parameters where our 
analysis is expected to hold qualitatively. 

In this work, we study the Josephson junction de- 













2 


scribed by a RCSJ model where the superconductors 
forming the junction either host Majorana modes at the 
interface or constitutes Dirac-like quasiparticles. In the 
former case, we show that in contrast to their counterpart 
in conventional junctions, the I-V characteristics display 
subharmonic odd Shapiro steps whose width vanishes for 
resistive junctions (in the limit where the junction ca¬ 
pacitance approaches zero)22122l \Ye provide an analyt¬ 
ical formula for the step width of both even and odd 
steps for such junctions, show that the analytic result 
matches exact numerics closely, and demonstrate, on ba¬ 
sis of this analytical result, that the behavior of this ra¬ 
tio is qualitatively different for Josephson junctions with 
and without Majorana bound states. In particular we 
demonstrate that the ratio of the odd and the even step 
width decrease exponentially with the junction capaci¬ 
tance for junctions with Majorana modes; in contrast, 
this ratio does not vary appreciably for conventional s- 
wave junctions. Thus it serves as a robust indicator for 
bound Majorana states in a JJ. Our analytical results, 
supported by numerical analysis, reproduces the phe¬ 
nomenon of absence of odd Shapiro steps in Josephson 
junctions with Majorana bound states as a special lim¬ 
iting property of resistive Josephson junctions; thus our 
work indicates that the absence of odd Shapiro steps, 
while sufficient, is not a necessary characteristics of Ma¬ 
jorana bound states in a Josephson junction. We also 
find that the I-V characteristics of junctions with Ma¬ 
jorana bound states show a qualitatively different devil 
staircase structure which is distinct from their s-wave 
counterparts. In particular, they display additional se¬ 
quence of steps which follow Farey’s sum rul^^ such se¬ 
quences are absent in I-V characteristics of conventional 
s-wave junctions. In the latter case, for junctions of su¬ 
perconductors hosting Dirac quasiparticles, we show that 
the width of the Shapiro steps display tt periodic oscilla¬ 
tory dependence on the barrier potential of the junction. 
Such a behavior is a direct consequence of the transmis¬ 
sion resonance phenomenon of the Dirac-like quasiparti¬ 
cles of the superconductors forming the junctions and is 
qualitatively distinct from conventional junctions hosting 
quasiparticles which obey Schrodinger equation. We note 
that our work shows that a RCSJ can act as phase sen¬ 
sitive detection device for both Majorana bound states 
in topological superconductors and Dirac like quasipar¬ 
ticles in a superconductor; it is therefore expected to be 
of interest to theorists and experimentalists working on 
both Majorana modes and Dirac materials. 


The plan of the rest of the work is as follows. In Sec. [Til 
we provide our analytical and numerical results for junc¬ 
tions which hosts subgap Majorana bound states. This 
is followed by Sec. [ml where we present our results on 
junctions which hosts Dirac-like quasiparticles. Finally, 
we summarize our main results, provide a discussion of 
experiments that can test our theory, and conclude in 
Sec. 113 



FIG. 1: (Color online) Schematic representation of the 
Josephson junction in the RCSJ circuit (see inset). The junc¬ 
tion has width L in the transverse direction and the barrier re¬ 
gion separating the two superconductors has a thickness d and 
is modeled by a potential Vq. Ij, Ic and Ir are the Joseph¬ 
son current, displacement current and quasiparticle currents 
respectively. 

II. JUNCTIONS WITH MAJORANA MODES 

The basic design of the circuit which we propose to 
serve as detector is shown in Fig. To analyze the 
property of this circuit, we first consider the Josephson 
junction. In our proposal, this comprises of two super¬ 
conductors with order parameters A/j (for x > d/2) and 
Ai (for X < —d/2) separated by a barrier region of width 
d {—d/2 < X < d/2) characterized by a barrier potential 
Vq as shown in Fig The superconductors can either 
be to polo gical superconductors with (effective) p-wave 
pairin^SlZI or superconductors with D irac-like q uasipar- 
ticles which has s-wave pair-potentiaP^E^^^lIIIl, jn this 
section, we analyze the former case in details. We note 
at the outset that the present analysis will hold for topo¬ 
logical superconductors in ID wire geometrjlSlH provided 
that the transverse dimension L is set to zero. 

Josephson junctions shown in Fig.[^ are known to sup¬ 
port localized subgap Andreev bound states which can be 
obtained as solution of the Bogoliubov-de Gennes (BdG) 
equation. For topological superconductors which support 
p-wave pairing, the BdG equations reads 

[{Hp + V(x))r^ + (Ap(x)T+ + h.c.)] i/p = Eipp (1) 

where P = R,L for the right and the left supercon¬ 
ductors, ipp = [ifjp'f-{x,k^^),ip'jj^{x,k^^)] is the two com¬ 
ponent BdG wavefunction, V{x) = Vq5{x) is the bar¬ 
rier potential (we take the limit of a thin barrier for 
which d ^ 0), Hp = hPk'^/{2m) — p denotes the dis¬ 
persion of the left and right superconductors with p be¬ 
ing the chemical potential, m the electron mass, and 

= -d"/. + fcjj- In what follows, we are going to as¬ 
sume Px-wave pairing and write A]^{x) = A^k^p/kp and 
Ap = AokxP/kp exjp{i(j)), where kp are the Fermi mo- 
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menta of the two superconductors, 4> is the phase dif¬ 
ference across the junction, Ag is the amplitude of the 
superconducting gap, and k^p is the x component of 
the Fermi momentum. The wavefunctions ipp satisfies 
the boundary condition ippix = 0) = 'iPr{x = 0) and 
dxi’Lix = 0 ) - dxtpRix = 0 ) = kFXi{k\\)ipL{x = 0 ), 
where Xii^w) = Xi = ‘^UQ/hvp{k^\) is the dimensionless 
barrier potential for a given transverse momentum of the 
quasiparticl es. T he localized subgap solutions of Eq. 
are given 

El = -Ao cos(</./2)/^1 + x2/4. (2) 

Note that for </> = tt, E = 0 and it has been shown that 
this state constitutes a realization of Majorana mode^. 
Since these subgap states are the only ones with (p depen¬ 
dent dispersion, one find the zero temperature Josephson 
current as 


where Pt — 1 requires the self-consistent solution of Eq. 
and Ij = dE 2 hyh/d(j). We have not attempted this in this 
work; however, we note that such a regime can always be 
obtained for long enough wire since 6 is exponentially 
suppressed in this regime leading to Pt ~ 1. An analysis 
of these conditions for resistive junctions can be found in 
Ref. dSl in the rest of this work, we shall assume Pt ~ 1 
and work with Eq. 

In what follows, we first obtain an approximate analyt¬ 
ical solution of Eq. I^in Sec. |II A| which demonstrates the 
existence of subharmonic odd Shapiro steps for junctions 
hosting Majorana bound states. This is followed by Sec. 

P where we carry out a detailed numerical study of Eq. 

lere the analytical results are verified and the devil 
staircase structure of the Shapiro steps is studied. 

A. Peturbative analytical Solution 


h{P)) = ‘^dEi/dP) = 


— kp 


dky eAg sin((/)/2) 

27r 2 ;iv^1 + x2/4- 


( 3 ) 


As noted in Ref.E^, the current is 47r periodic and a sub¬ 
stitution (j) —>■ 2eVt/h in the presence of a bi as voltage V 
leads to fractional AC Josephson effectP^^^fl, 

We now use Eqs. [^to obtain the response of an RCSJ 
circuit constructed out of the superconducting junctions 
discussed above. The RCSJ model, shown in the inset 
of Eig. include a resistive component to take into ac¬ 
count dissipative process into account which may occur, 
for example, due to quasiparticle tunneling and a shunt¬ 
ing capacitance C which takes into account the displace¬ 
ment currents due to possible charge accumulation in the 
lead^^. The current phase-relationship for this model, in 
the presence of an external radiation, is given by^SM! 


^ + + = 7//, + Asin(wt)//e, (4) 


In this section we consider perturbative analysis of Eq. 
for an unconventional Josephson junction comprised 
of superconductors hosting subgap Majorana states with 
E = Ej cos{(j>/2). We begin from the equation of such 
a junction given by Eq. and analyze this equation for 
uj,l3u),A » 1. The key point regarding this analysis 
is the observation in the regime mentioned above it is 
possible to expend (p asP 

OO 

</) = I = (5) 

n n=0 

where /g is the applied current and e <C 1, and for 
n > 0 determined self-consistently from the condition of 
absence of additional dc voltage: lim-r-^oo (pn = (P- 
The equations for (pn can be obtained by equating 
terms in the same order of e. The procedure is standard 
and yields 


where Ic is the critical current of the junction, A and 
Lo are the amplitude and frequency of the external ra¬ 
diation, Ij{(p) = h for superconducting junctions with 
Majorana fermions, (3 = ^hji^elcEG-Ct^, R and Cq de¬ 
note the resistance and capacitance of the junction, and 
we have scaled t — >■ t/r, where t = ^yhCo/{ 2 eIc)- 

Before proceeding further, we note that Eq.[^ holds in 
the ideal limit where the two subgap Majorana modes 
at the two ends of each wire do not interaciP. Such an 
interaction leads to hybridization of amplitude 6 between 
the two Majorana branches: Eihyb = ±\/J^ + E^^. 
The hybridization amplitude <5 arising from such inter¬ 
action is exponentially suppressed for long wires S ~ 
exp[—E/(2^)], where the coherence length for ID 
nanowires, depends on the spin-orbit coupling strength 
of the wire. In the presence of an external voltage, this 
leads to a Landau-Zenner tunneling probability Pt = 
exp[—27rJ^/(Ejfi(^)], where Ej ~ Ag is the maximal 
Josephson energy. The manifestation of Jtt periodicity is 
evident when ~ 1. This occurs when 6'^/{Ejh(p) <C 1; 
a complete determination of frequency and voltage range 


(pn + = fn{t) +4, /o = As\a{lot), 

fi = -sin[4/2], f 2 = (piCOs[(po/2]/2 (6) 

Note that the n = 0 equation represents the autonomous 
I-V curve of the junction and is independent of the non¬ 
linear sinusoidal term. 

To solve these equations, we note that this represent a 
linear first order differential equation in (pn] consequently 
we define = (pn and write 

yn+ Pyn= In + fnit). (7) 

These equations admit a solution 

yn(t) = In + e~^* [ fn{t')dt', 

Jo 

4>n{t) = [ yn{t')dt' (l)n{0). (8) 

Jo 

For n = 0, this yields 

A 

(j)o{t) = (j)' iQt-\ - sin{ujt + ao), (9) 

CU'J 


















4 


where ao = arccos(a;/ 7 ), 7 = \//3^ + and </>' is the 
DC phase of the junction. The supercurrent at this order 
is given by 

4°) - sin(<^o(t)/2) = (10) 

00 

= Im ^ j^( 3 ,)g*([/o/ 2 +na;]i+nao+^'/ 2 )^ 

n——oo 

where x = A/{2'yL>j). Thus the Shapiro steps occur 
when the AC component of the supercurrent vanish: 
/o = 2|n|a; (even steps). The width of the even step 
can be read off from Eq. W^even = = 2J„(a:). 

These yield expressions for the width of harmonic steps. 
Note that in contrast to the conventional junctions where 
la ~ sin[(/)o], only even harmonic steps occur for junctions 
which host subgap Majorana steps. 

Next, we obtain the solution for (j)i. Substituting Eq. 
11 in Eq. we find after some straightforward algebra 


(/)i = ^ J„(x)( 7 „w„) ^ cos{u!nt + nao + So + ncj)'/2). 


n— — oo 


( 12 ) 

where a;„ = Io/2 + nw, Sn = arccos(a;„/ 7 „), and 7 „ = 
y/oJn + At this order, the supercurrent is given by 

4^^ ^ ^</>iWcos(0o(4/2) 


^ ^ All (x)(Ai 2 (^) (27^1 ) 

ni ,712——cxD 

X sm{uJnA + ni{ao + 4>'/2) + 

X COs{uJ„^t + 712 (ao + ((> 72 ) + SnJ 
7 ^ (x) t/yi2 (^) (47ni ) 


-1 


sin([w„i +UJn2]t + [ni + 7i2](ao + 7/2) + (^m) 


+ sm([a;„i - Wnjt + [ni - 772 ](ao + 7/2) + <Jni) 

At this order, we find that there are additional steps 
in the DC component of the supercurrent; these steps 
occur at | 77 i + 772 |w = Iq for which the first of the two 


terms in the right side of Eq. 14 become independent of 
time. A set of these steps occur at {ni + 772 ) = 277 —1 
for integers n = 1,2... and constitute the odd Shapiro 
steps. Thus we find that the odd steps for a junction 
of superconductors hosting Majorana ground states are 
necessarily of subharmonic nature. The width of these 
steps can be read off from Eq. to be 


VEodd = A/° 


dd 

sn 


E 


dni (^)d2n—1 —ni (^) 

2[({277 — 1 — 77i}a;)^ + 


<14) 


We note that when Cq —>■ 0, /3 —>■ 00 and the subharmonic 
steps vanish leading to the result that only even har¬ 
monic Shapiro steps exists for resistive Josephson junc¬ 
tions hosting subgap Majorana steps. Thus our analysis 


reproduces the absence of odd Shapiro steps in Joseph¬ 
son junc^ns with Majorana bound states as a special 
cas^iSlMM]^ We also note that these odd steps have a 
completely different origin than the analogous steps dis¬ 
cussed in Ref. [53] since they occur without any ~ sin(^) 
dependence of Ij. Einally, we note that the ratio of the 
77 *'' even and the adjacent odd Shapiro steps for these 
junctions are given by 




= 


2J„(x) 


W^odd 


2-^ni 2[({2n—1—ni}cj)^+/3^] 


(15) 


(13) 


In the next section, we shall compare the analytical ex¬ 
pression (Eq. 151 with numerical results obtained by ex¬ 
act numerical solution of Eq. 


B. Exact numerical results 

To compute the I-V characteristics, we study the tem¬ 
poral dependence of E = h(j)/{2e) obtained by numerical 
solution of Eq. I^as a function of time for a fixed bias cur¬ 
rent I. The dc co mpone nt of the voltage is obtained by 
standard procedur d^^l'"' ! from V and plotted as a function 
of I to generate the I-V characteristics. 

The central results that we obtain from this analy¬ 
sis are as follows. First, for topological superconductors 
hosting Majorana subgap states, we find that for a signif¬ 
icant range of external coupling and in the underdamped 
region /3 < 1, show both even and odd Shapiro steps 
as expected from the analytical results obtained in Sec. 
II A[ The even steps at E = 2ntiuj/e are enhanced com¬ 


pared to their odd counterparts at E = {2n + l)tujjje as 
shown in Figs. and for underdamped (/3 < 1 ) and 
overdamped (/3 > 1) regions respectively. Fig. shows 
an analogous plot for the s-wave superconductors. The 
dominance of the even steps over the odd ones is charac¬ 
terized by 7^1 = ?7 (Eq. 15). A plot of 7 y as a function of /3 
shown in Fig. demonstrates that 77 ^ exp(0.3/3^) for 
junctions with Majorana modes. We also note that the 
theoretical result for 77 obtained from Eq. provides a 
near-perfect match with the exact numerics demonstrat¬ 
ing the accuracy of the analytical solution over a wide 
range of /?. We further note that the behavior of 77 as a 
function of /? is in complete contrast to its counterpart for 
s-wave superconductors where 77 does not vary apprecia¬ 
bly with /3 as shown in Fig. W- Thus the exponential de¬ 
pendence of 77 on the junction capacitance C constitutes 
a phase sensitive signature of the presence of the Majo¬ 
rana modes. Note that it is generally expected that only 
even Shapiro steps occur in I-V characteristics of Joseph¬ 
son junctions which supports Majo rana mo des due to Jtt 
periodicity of the Josephson curren-rfi^ISlIlSl jg ^ con¬ 
sequence of analysis of the problem in the limit of zero 
junction capacitance C —>■ (E^EZEll jj ^ leading 

to vanishing of the odd steps. However, our study show 
that underdamped RCSJ of such unconventional Joseph¬ 
son junctions with finite C can display both even and 
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FIG. 2: (Color online) Top panels: CVC of p-wave Josephson 
junction. Fig[^ shows the CVC in the underdamped region 
(/3 = 0.2) and for A = 20 and lu = 2. while Fig. |2> represents 
the overdamped region (/3 = 1.2) with A = 5 and w = 1 (in 
appropriate dimensionless units; see text). Bottom Panels: 
Fig.§ shows the CVC of s-wave junction (all parameters are 
same as in Fig. Fig.§i shows the ratio of the widths of 
the even and odd Shapiro steps, rj = II4ven(2a;)/Wodd(a)) for 
the p-wave and the s-wave junctions (inset) as a function of 
f] for j 4 = 10 and to = 3. For p-wave, rj ~ gO.sifi while for 
s-wave, rj does not vary with j3. Figs. [^a), (b), and (c) have 
V and I scaled in units oih/ (re) and Ic respectively. The red 
solid (blue dashed) curves in panels (a) and (c) correspond 
to data for increasing (decreasing) current sweeps; these data 
coincides in the overdamped region as shown in panel (b). 


odd steps. Thus the presence of odd Shapiro steps do 
not necessarily signify the absence of Majorana modes 
specially if the RCSJ is underdamped. 

Another qualitative difference between the I-V charac¬ 
teristics of Josephson junctions hosting Majorana states 
with their conventional counterparts occur in the devil 
staircase structures of the Shapiro steps occurring within 
a fixed bias current intervals in these junctions. Such 
steps are known to occur for conventional s-wave junc¬ 
tions and their voltage-frequency relation can be repre¬ 
sented by a continued fraction as 

V = [Af ± l/(n± (1 /to ± l/(p± 1...)))] w, (16) 

where N,n,m,p are integer^^. The fractions involving 
N are termed as first level fractions, those with N and n 
are second level fractions, and so on. To distinguish be¬ 
tween s-wave and p—wave junctions, we investigate the 
structure of these steps for both p— and s—wave junc¬ 
tions as shown in Figs. andWe find that for s-wave, 
only second level subharmonics are present and V = Olu 
is approached from below with steps occurring ai V = 
N{l—l/n)uj for different N and n. In contrast, for p-wave 
junctions, one finds that in addition to the steps which 
occur for the s—wave junctions, there are additional steps 



0.575 0.58 0.585 |0.59 0.595 


FIG. 3: (Color online) Plots of the self-similar structure for 
p-wave (A = 0.6) s-wave {A = 0.8) Josephson junctions for 
D = 0.7, and ui = 0.5. These additional fractions marked with 
arrows pointing to the right belong to the additional sequence 
characteristics of Josephson junctions with Majorana snbgap 
states and obey Farey sum rule. V and I are scaled in units 
of hj{eT) and Ic respectively; see text for details. 


corresponding to V = [(Af =F 2)][1-I-l/(n± 1 /to)]u; (Figs. 
[^and|^. This difference in structure lead us to hypothe¬ 
size that in contrast to s—wave Josephson junctions, the 
continued fractions for p—wave shows additional fractions 
all of which are consistent with the Farey sum rule; z.e., 
the widest phase-locked region (the widest step) between 
any two resonances p/q and m/n for these fractions is 
given by {pm) / [q + n). We therefore suggest that the 
presence of these additional specific continued fractions 
may be considered as a signature of Majorana modes. 

Two specific examples of this phenomenon is presented 
in Fig. and In Fig. we find that the continued 
fraction V = {N — l/n)a; with N = 6 which appears 
in s—wave Josephson junctions at /? = 0.2, ui = 0.5 
and A = 0.8, is also appeared in the p—wave junction 
at smaller amplitude of radiation A = 0.6. However, 
as clearly demonstrated in Fig. the steps for the p— 
wave junction constitutes an additional sequence of con¬ 
tinued fraction V = [{N — 2)][1 + l/(n -|- l/m)\uj with 
N = 6 , n = 2 and for several to . Further, as shown in 
Ref. Uni the continued fraction V = {N -|- l/n)uj with 
N = 6 which appears in s—wave Josephson junctions at 
A = 0.9. In contrast, for p—wave Josephson junctions, as 
shown in Fig. there is an additional continued fraction 
V = [(A^-|-2)/2][l-|-l/(n—1 /to)]w with N = 6 and n = 2. 
In both cases, additional steps occur when Majorana sub¬ 
gap states are present; thus, although we have not been 
able to find a precise analytical expression relating the Jtt 
periodicity of Ij{(j)) and these additional steps, we con¬ 
jecture that these additional steps are a consequence of 
the subgap Majorana bound states of the superconduc¬ 
tors forming the junctions. 
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FIG. 4: (Color online) Same as in Fig. but with different 
amplitudes of external radiation [ p-wave (A — 0.77) s-wave 
{A = 0.9)]. Note that the additional sequence of steps for p— 
wave persists. 


III. JUNCTIONS WITH DIRAC FERMIONS 

For superconducting ju nction s hosting Dirac quasipar¬ 
ticles such as in graphen^^^Ull^ pair-potential has s- 
wave symmetry: Ap = Aq and A^ = Aoexp(z(/)). Here 
pairing occurs between electrons with opposite spin and 
momenta; in graphene this necessitates pairing between 
electrons of K and K' valleyJi^. This the BdG equa¬ 
tions are described in terms of 4 component wavefunc- 
tions ijj = V’s/) and is given by 

[{H' - V{x))t 3 -h (A(a:)r+ -|- h.c.)]-)/^ = E-ijj, (17) 

where A, B denotes sublattice indices = 

hvp{—iaxdx ± (Xyky) — n, t and a denote Pauli matrices 
in valley and pseudospin (sublattice) spaces respectively, 
pL is the Fermi energy, vp is the Fermi velocity of the 
Dirac quasiparticles described by H' and -!-(—) sign cor¬ 
responds to electrons in K{K') valley. The pair-potential 
takes the form A (a;) = AL9{d/2 + x) + Ap6{x — d/2) and 
V{x) = Vo9{d/2 — x)9{d/2 + x). The localized subgap 
Andreev states is obtained by demanding the continuity 
of '0 at a; = Ad 12 and is given, in the thin barrier limit 
(Vb oo, d —>■ 0, aTidVod/hyp = X 2 ) 


19 shows that the 


where Iq = eAQEpL/{2H^7rvp). Eq. 

Josephson current is an oscillatory function of the di¬ 
mensionless barrier strength for such junctions. We note 
that whereas the form of Eqs. and are generic for 
s—wave conventional tunnel junctions, the oscillatory de¬ 
pendence of T on X is a consequence of Dirac nature of 
graphene quasiparticles and is not observed in junctions 
made of conventional superconductor^^. 

To chart out the I-V characteristics of the supercon¬ 
ducting junctions which host such Dirac quasiparticles, 
we analyze Eq. numerically with Ij = I 2 and obtain 
the corresponding Shapiro step structure. The proce¬ 
dure followed here is identical to the one charted out in 
Sec. |IIB[ We find that the Shapiro step structure in the 
I-V characteristics is same as the conventional s-wave su¬ 
perconductor displaying harmonic odd and even steps as 
shown in Fig. [^. However, the width of these steps, W, 
vary with the dimensionless barrier potential X 2 in an os¬ 
cillatory manner as shown in Fig.[^. This is in complete 
contrast with the dependence of W in the conventional 
junctions where the step widths are monotonically de¬ 
creasing function of the barrier potential. This behavior 
of W can be qualitatively understood as follows. In a 
RCSJ, W, can be related to the magnitude of Ij which, 
in turn, depends on the transparency of the junction: 
W {1 + for conventional junctions. In con¬ 

ventional junction, increase of the barrier potential xi 
leads to a monotonic decrease of the transparency; con¬ 
sequently, W decrease monotonically with increasing xi- 
However, for a RCSJ made out of Dirac materials, the 
transparency of the junction, T{ky, X2), is an ir periodic 
oscillatory function of the dimensionless barrier strength 
X 2 with maxima at x = due to the transmission res¬ 
onance condition of the Dirac quasiparticleJi^. Conse¬ 
quently, one expects Ij and hence W to oscillate with 
X 2 . This expectation is corroborated in Fig. [^, where 
the TT periodic oscillation of the step width is plotted as a 
function of X 2 - We note that such a oscillatory behavior 
is a direct manifestation of the transmission resonance 
condition of the Dirac quasiparticles; it thus provides a 
qualitative distinction between Josephson junction host¬ 
ing Schrodinger and Dirac quasiparticles. 


IV. DISCUSSION 


E 2 = ±A3^Jl-T{ky,X2)Bin\cl>/2) 


(18) 


where T{ky,X 2 ) = cos^( 7 )/(l — cos^(7) sin^(x 2 )) is a 
measure of transparency of the junction and sin( 7 ) = 
hvpky/pL. Note that Eq. 18 in contrast to Eq. 

is 2p periodic in </). The corresponding Joseph¬ 
son current at zero temperature is given by l 2 { 4 ') = 
V ^’^^i'y)^E 2 /d(l) and leads to 

. 77/2 

I2 = IqAo / d7COs(7)sin((;j)r(7, X2)/|£'2| (19) 

J —7^12 


In this work we have studied the I-V characteris¬ 
tics of a RCSJ where the superconductors making up 
the junction either hosts subgap Majorana bound states 
or have Dirac-like character of the Bogoliubov quasi¬ 
particles. The former set of junction s o ccur for p- 
wave superconductorJj^ or ID nanowire^^ with strong 
spin-orbit coupling and transverse magnetic field while 
graphene superconduction junctions provide an example 
of the latter class. We find that the I-V characteristics 
of RCSJs for each of these classes of junctions are qual¬ 
itatively different from their conventional counterparts. 
Thus such junctions may serve as phase sensitive detec- 
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FIG. 5: (Color online) (a) CVC of graphene Josephson jnnc- 
tion with A = 1 and w = 4 for several values of x indicating 
the variation in the width of the main Shapiro step at uj. (b) 
Plot of the width of the Shapiro step at a; as a function of x 
showing TT periodic oscillatory behavior. 


tors of Majorana and Dirac fermions realized using su¬ 
perconducting platforms. 

For junctions hosting subgap Majorana states, we find 
two essential characteristics which are qualitatively dif¬ 
ferent from their s-wave counterparts. First, the odd 
Shapiro steps are subharmonic in nature; the ratio of 
their width with that of adjacent even Shapiro is a de¬ 
creasing function of the junction capacitance. This is in 
contrast to the conventional s—wave junctions where the 
ratio is largely independent of C. We note that our re¬ 
sult in this regard shows that the absence of odd Shapiro 
steps is a sufficient condition for having subgap Majorana 
modes; however it is not necessary since such an absence 
requires, in addition to the presence of the Majorana 
modes, resistive Josephson junctions. Our result thus 
constitute a generalization of detection criteria for Ma¬ 
jorana modes realized using superconducting platform. 
Second, we find that the devil staircase structure of the 
Shapiro steps in Josephson junctions with Majorana sub¬ 
gap states involves additional sequences which satisfies 
Farey sum rule. This feature, as shown in our work, con¬ 
stitutes a qualitative difference between Josephson junc¬ 
tions with and without Majorana subgap states. 

For junctions with Dirac quasiparticles, we find that 
even with s-wave symmetry, the Shapiro step width is 
a TT periodic oscillatory function of the barrier potential 
of the junction. We trace the origin of this phenomenon 
to the transmission resonance of the Dirac-Bogoliubov 
quasiparticles in such superconductors and demonstrates 
that the oscillatory behavior is a qualitatively distinct 
signature of Dirac nature of the superconducting quasi¬ 
particles. 

The numerical estimate of typical frequencies at which 
the devil staircase structure can be obtained as follows. 
For standard experiments Ic ~ InA and Co — IpF. Using 
these numbers one can estimate, ujp = •\/2eJc/(^C q) — 


IGHz. In all the plots, that we have used uj ranges 
between 0.5wp ~ 0.5GHz to 2ujp ~ 2GHz. In particu¬ 
lar, the devil staircase structure is seen at an external 
radiation frequency of 0.5GHz. The self-similar struc¬ 
ture is seen at energy range of 5 — Gw which is around 
2 — 3GHz. In this context, we note that the required 
frequency range is large small enough to avoid possi¬ 
ble smearing due to 27r periodicity arising due to mul¬ 
timode effect^^. In addition, this estimate holds for 
zero-temperature analysis; however it is expected to be 
qualitatively accurate for kgT Ag when quasiparti¬ 
cle poisoning and thermal decoherence rates do not play 
a significant role. Also, we note that it is possible to 
model the dissipation and noise in the junction by as¬ 
suming it to be coupled through a thermal bath using the 
standard Caldeira-Leggett formalisnP^, the Langevin or 
saddle point equation corresponding to that analysis at 
low temperature, where the effects of white noise can be 
ignored, reduces to Eq. of our work with the resistive 
term being renormalized by the coefficient of dissipation. 
This formalism also allows for study of effect of quantum 
fluctuations and noise in such junctions beyond saddle 
point approximation which is left as a possible subject of 
future study. 

The experiments to test our theory involves on mea¬ 
surement on RGSJ under an applied radiation with defi¬ 
nite amplitude A and frequency w. Such experiments are 
rather standard for s—wave junction^^ more recently, 
such experiments have been performed for ID Majorana 
nanowire setup with resistive junctionJi^. Our specific 
suggestion involves measurement of ?7 as a function of 
the effective junction capacitance for JJS with Majorana 
bound states in a circuit with finite capacitance which 
can be modeled by a RGSJ; we predict the presence of 
subharmonic odd Shapiro steps for such junctions whose 
width depends on the junction capacitance and lead to 
an exponential dependence of r] with the junction capaci¬ 
tance (See Fig. W)- In addition, we suggest the presence 
of an additional sequence in the devil staircase structure 
of the Shapiro steps. For junctions with Dirac quasipar¬ 
ticles, which can be made with graphen^^, we predict 
that the width of the Shapiro steps will display tt pe¬ 
riodic oscillatory dependence with the junction barrier 
potential. 

In conclusion, we have studied RGSJ Josephson junc¬ 
tion circuits and have shown that they can serve as phase 
sensitive detectors for both Majorana and Dirac quasi¬ 
particles in such junctions. We have charted out the 
properties of such junctions which are qualitatively dis¬ 
tinct from their s-wave counterparts and have suggested 
experiments which can test our theory. 
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